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ABSTRACT 

Aims. We characterize the molecular gas content (ISM cold phase) using CO emission of a redshift-limited subsample of isolated 
galaxies from the AMIGA (Analysis of the interstellar Medium of Isolated GAlaxies) project in order to provide a comparison sample 
for studies of galaxies in different environments. 

Methods. We present the l2 CO(l-0) data for 273 AMIGA galaxies, most of them (n = 186) from our own observations with the IRAM 
30m and the FCRAO 14m telescopes and the rest from the literature. We constructed a redshift-limited sample containing galaxies 
with 1500 km s~ l < v <5000 km s -1 and excluded objects with morphological evidence of possible interaction. This sample (n = 173) 
is the basis for our statistical analysis. It contains galaxies with molecular gas masses, Mu 2 , in the range of ~ 10 s - 10'° M Q . It is 
dominated, both in absolute number and in detection rate, by spiral galaxies of type T = 3 - 5 (Sb-Sc). Most galaxies were observed 
with a single pointing towards their centers. Therefore, we performed an extrapolation to the total molecular gas mass expected in the 
entire disk based on the assumption of an exponential distribution. We then studied the relationships between M H , and other galactic 
properties (Lb, D 2 5 , Lk, Lfir, and Mm) • 

Results. We find correlations between M Hl and L B , D^ 5 , L K , and L FiR . The tightest correlation of M Hl holds with L F i R and, for 
T = 3 - 5, with Lk, and the poorest with D 2 5 .The correlations with Lfir and Lk are very close to linearity. The correlation with L B 
is nonlinear so that M H ,/L B increases with L B . The molecular and the atomic gas masses of our sample show no strong correlation. 
We find a low mean value, log(MH 2 /Mm) = -0.7 (for T = 3 - 5), and a strong decrease in this ratio with morphological type. The 
molecular gas column density and the surface density of the star formation rate (the Kennicutt-Schmidt law) show a tight correlation 
with a rough unity slope. We compare the relations of Mh 2 with L B and Lk found for AMIGA galaxies to samples of interacting 
galaxies from the literature and find an indication for an enhancement of the molecular gas in interacting galaxies of up to 0.2-0.3 
dex. 

Key words, galaxies: evolution - galaxies: interactions - galaxies: ISM - radio lines: ISM - radio lines: galaxies - surveys 



• rH ! 1- Introduction vironme nt affects the amount of the molecular gas. Galaxies in 

■ clusters dKennev & Yound[l98 9; Bosel li et al.ll 19971 Scott et al. 

' ma J° r and pasting debate in astronomy involves the in prep .) and groups fcdes-Monteneg^Saffi998; LeoneUD 

Ctf ; relative roles of nature and nurture in galaxy forrmtion 1 1998) seem tQ haye a normal molecular gas contentj even though 

and evolution (e.g. Sulentic 1976; j Larson & Tinslev | [1978j they can be highl y deficient in a tomic gas . On the other hand, 

Uoseph & Wnght| | l985HBushouse| | l987|) . Although it is broadly SQme authors , Br aine & Combes 1993; Combes et al. 1994; 

accepted that galaxy evolution strongly depends on the envi- Casasola et al 2004) find an enhanced mo i e cular gas content 

ronment, the quantitative effect of "nurture" on certain galactic in interacti ga i axies , in contrast to the results of iPerea et all 

properties is still a matter of debate. ([1993), who concluded that the molecular gas content is not af- 

The molecular gas content is an important quantity of a fec ted by interaction in strongly interacting pairs or Virgo cluster 

galaxy because it is directly related to its capacity for star for- galaxies 

mation (SF). We still need to determine, however, how the en- „ , .„ , , , , , £ , 

To clarity the role played by the environment, a well-defined 

* t Full Tables 1, 4 and 5 are only available sample of isolated galaxies is needed to serve as a zero level 

in electronic form at the CDS via anonymous ftp for studles dealing with the effect of interactions. Most pre- 

to cdsarc.u- s trasbg.fr (138.7 9.128.5) or via vious studies investigating the p roperties of molecular gas in 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J /A+A/ and isolated and interacting galaxies ( Solomon & Sage 1988; Sage 

from |http : //amiga . iaa . esTT ll993blBoselli et alJl997tlNishivam a & Na kail200lLlHelfer et al. ' 
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12003b iLerov et al.1 120091) h ave generally not defined any very 
clear criterion for isolation. iPerea et aD d 1 9971) carried out a CO 
study comparing isolated and interacting galaxies. Their sam- 
ple of isolated galaxies is composed of 68 galaxies from var- 
ious sources, selected in a much less rigorous way than the 
present study and biased towards infrared-luminous objects. The 
only survey explicitly focusing on isolated galaxies, and in par- 
ticular on galaxies from the C atalogue of Isolated Galaxies, is 
the one by ISautv et al] (120031) . They present the CO data of 99 
optically-selected spiral galaxies with recession velocities up to 
14000 km s and briefly compare the properties of the molecu- 
lar gas mass to the blue luminosity and atomic gas mass. A de- 
tailed analysis of the properties of that sample is, however, not 
presented there. The larges t previous CO sur vey was the FCRAO 
Extragalactic CO survey (lYoung et al.l 1 1995b observing ~ 300 
nearby galaxies. The major difference with respect to the present 
study is that it did not consider the isolation of the galaxies as a 
criterion. Furthermore, it only contained bright galaxies (either 
m B ,coiT < 13 mag, or F 60fm > 5 Jy or Fioo^ > 100 Jy) whereas 
our samples also includes fainter objects. 

The project AMIGA ('Analysis of the interste l lar Me dium 
of Isolated GAlaxies", IVerdes-Montenegro et al.1 [2005 ) was 
started to provide such a reference sample by characteriz- 
ing the properties of the interstellar medium (ISM) and star 
formation (SF) in isolated galaxies in the local Universe. 
It is based on the Catalogue of Isolated Galaxies (CIG 
iKarachentseval [l973h which is composed of 1050 galaxies lo- 
cated in the Northern hemisphere. The A MIGA project is 
presented in IVerdes-Montenegro et al.l d2005l) . A considerable 
amount of work has been done since then in order to refine 
the sa mple. This work includes the re vision of all CIG posi- 
tions dLeon & Verdes-Monten egro 2003), the determination of 
POSS2-based morphologies and the i dentification of gal axies 
showing signs of possible interaction (Sule ntic et al.ll2006h and 
the reevaluation and quantification of the degree of isolation 
dVerlev et al.ll2007allbl) . The results of this project consistently 
find that the AMIGA galaxies have the lowest SF activity as 
well as the lowest presence of Active Galactic Nuclei (AGN) 
in the local Universe. This is obtained bot h from the far-infrare d 
(FIR) luminosity derived from IRAS data (Lisenfeld et al. 2007), 
and the radio continumm emission dLeon et al.l 120081) . which 
are both SF tracers. The rate of AGN candidates, derived from 
IRAS colors and radio continuum emission of the AMIGA 
galax ies is lowest compa red to similar studies from the litera- 
ture (Sab ater et al.l [2008). Optical photometric analysis of Sb- 
Sc galaxies in the AMIGA sample showed that most galaxies 
have pseudo-bulges instead of classical bulges, and a compar- 
ison with samples of spiral galaxies selected without isolation 
criteria revealed that the isolated galaxies tend to host larger 
bars, are more symm etric, less concentrated and less clumpy 
dDurbala et al.ll2008l) . These findings strongly support that the 
AMIGA sample represents the most isolated galaxies in the lo- 
cal Universe where secular evolution i s dominant. Espad a et al. 
(in prep.) study the HI content and dEspada et al.l |20 1 lb found 
the smallest fraction of asymmetric HI profiles in the AMIGA 
sample when compared with any sample yet studied. 

The revised AMIGA sample is re asonably complete (~ 80- 
95%) down to niBcorr ^ 15.0 mag (IVerdes-Montenegro et alJ 
120051) and it is currently one of the largest sample of nearby iso- 
lated galaxies in the Northern hemisphere. It consists of galaxies 
whose structure and evolution have been driven largely or en- 
tirely by in ternal rather than by external forc es at least during the 
last 3 Gyr dVerdes-M ontene gro et al.ll2005l) . The data are being 
released and periodically updated at http:// amiga.iaa.es where 



a Virtual Observatory interface with different query modes has 
been implemented. 

In the present paper we present and analyze CO observations 
of a redshift-limited subsample of this catalogue. The goal is to 
characterize the properties of the molecular gas, traced by CO, 
of isolated galaxies and to provide a reference sample for studies 
investigating the role of the environment. 

2. The sample 

For a study of the molecular gas content we had to restrict 
the number of galaxies, since observation of the entire opti- 
cally complete sample (n ~ 700) required too much telescope 
time. We chose to build a redshift-limited subsample by select- 
ing galaxies with recession velocities in the range of 1500-5000 
km s _1 . The completeness limit of of the AMIGA sample of 15 
mag corresponds to blue luminosity of logl0(LB/L o ) = 8.55 and 
logl0(LB/L o ) = 9.60 at the distances derived for these veloc- 
ities with a Hubble constant of 75 km s _1 Mpc -1 . The range 
was chosen in order to avoid (i) very n earby galaxies for w hich 
the condition of isolation is not reliable dVerlev et al.ll2007ah and 
(ii) distant galaxies which are difficult to detect in CO. The re- 
striction in velocity provides us with a sample probing a defined 
volume in space. 

There are 278 galaxies in this velocity range in the CIG. We 
have CO data for 201 of these objects, mostly from our own ob- 
servations (180 galaxies) with the 30m telescope of the Instituto 
de Radioastronorma Milimetrica (IRAM) at the Pico Veleta and 
with the 14m Five College Radio Astronomical Observatory 
(FCRAO), and the rest from the literature. We then excluded 
those galaxies that were identified by us in a visual inspection 
of optical images as having signs of a possible present or past 
interaction (see description of Table Q] for more details on the 
criteria). This leaves us with 173 isolated galaxies with CO data 
in the velocity range between 1500 and 5000 km s . We refer 
to this sample as the redshift-limited CO sample and we will use 
it for the statistical analysis throughout this paper. 

Additionally, we have CO data for 72 galaxies outside this 
velocity range. Six galaxies are from our own observations (with 
velocities between 5000 and 5500 km s _1 ) and the rest is from 
the literature. Thus, in total, we have 12 CO(1-0) data for 273 CIG 
galaxies. We refer to this sample as the total CO sample, and list 
the corresponding data in Tables[T]and|5] but we do not use it for 
any statistical analysis. 

Recently, an update of the basic properties of the galaxies 
in the AMIGA sample was carried out for the blue magnitude, 
optical isophotal diameter D25, velocity, and morphology and in- 
teraction degree based on higher resolution images (from SDSS 
or our own images). The details are described in Espada et al. 
(in prep.). In order to provide a self-contained data set for the 
present paper, we list in Table Q] the basic data relevant for the 
total CO sample. The columns are: 

1. CIG: Entry number in the Catalogue of Isolated Galaxies 
(CIG). 

2. Dist: Distance in Mpc, based on Ho = 75 km s 1 Mpc . 

3. Vel: Recession velocity in kms -1 . 

4. D25: Isophotal optical major diameter at the isophotal level 
25 mag arcsec -1 in the B-band, D25, in arcmin from the Lyon 
Extragalactic Database (LEDA). 

5. i: Inclination from LEDA. 

6. T(RC3): Morphological type T(RC3) as determined by our 
morphological revision, giv en in the RC3 numerical scale 
dde Vaucouleurs et aD 19911) . 
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CIG 


Dist 


Vel 


D„ 


i 


T(RC3) 


Inter 


log(L B ) 


log(L F iR) 


log(L K ) 




[Mpc] 


[kms-»] 


['] 


[°] 






[Lq] 


[L Q ] 


[£k,o] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


96 


7299 


1.39 


65 


5 


1 


10.57 


10.28 


11.23 


4 


31 


2310 


3.29 


84 


5 





10.36 


10.08 


10.98 


6 


61 


4528 


0.70 


65 


7 


1 


9.80 


9.82 


10.21 


10 


63 


4613 


1.05 


63 


5 





9.78 


<9.44 


10.09 



The full table is available in electronic form at the CDS and from http://amiga.iaa.es 



7. Inter: Code for morphological signs of a possible interaction. 
0: No signs of large distortions which cannot be clearly ex- 
plained by e.g. dust. 

1 : At least one of the following signs of possible interaction: 
asymmetric, lopsided, distorted, co nfirmed pair in the Nasa 
Extragalactic Database (NED) or in lVerlev etall d2007bl) . in- 
tegral sign shape, warp, tidal features (tail, bridge, shell). 
2: Merger-like morphology or superposition of two galaxies. 

8. log(LB): Decimal logarithm of the blue luminosity, Lb, in 
units of solar bolometric luminosities. We calculated Lb, 
from the flux, fs, as L^= vfs, where v is the central fre- 
quency of the blue band. 

9. log(LFiR): Decimal logarithm of the far-infrared luminos- 
ity in units of solar bolometric luminosities, taken from 
iLisenfeld et alJ (120071) and adapted to the revised distances 
used in the present paper. Lfir is is computed from the IRAS 
fluxes at 60 and 100 yum, F6o and Fioo), as log(LpiR/Lo) = 
log(FIR) + 21og(D) + 19.495, where D is dist ance in Mpc 
and F IR = 1.26 x 10" 14 (2.58F 60 + Fioo) Wm" 2 dHelou et alJ 
Il988h . 

10. log(LK): Decimal logarithm of the luminosity in the K-band, 
in units of the solar luminosity in the A^-band (Lk,o = 
5.0735 x 10 32 erg s~'), calculated from the extrapolated mag- 
nitude in the K$ (2 .17 um) band fro m the 2MASS Extended 
Source Catalogue (Jarrett et al. 2000). The magnitudes were 
available for 250 galaxies of our sample. We calculated the 
Ks luminosity, Lk, from the total (extrapolated) K$ flux, /k, 
as Lk = K7k(v), where v is the central frequency of the Ri- 
band. Lk is a good measure of the total stellar mass. 



In FigJTJwe present some of the basic characteristics of the 
CO sampl es and compare them to those of the optically complete 
sample (I Verde s-Mon tene gro et al. 2005). The latter (n - 712) is 
composed of CIG galaxies with wb in the range of 1 1 - 15 mag 
and is 80 - 95 % complete. The main difference between the CO 
and the optically complete sample is the larger spread in velocity 
of the latter (FigQJ). The CO samples are at a lower velocity, 
especially the redshift-restricted sample. This leads to a slightly 
lower optical luminosity (FigJTJ?) of the CO samples since the 
number of luminous objects is higher at larger distances due to 
the Malmquist bias. 

The distribution of the morphological types (FigJTJ;) is very 
similar in all the three samples. All samples are dominated by 
spiral galaxies. The relative number of early type galaxies (T = 
(— 5)-0) is 8% for the optically complete sample and 13% for the 
CO redshift-restricted sample. The distribution is peaked around 
galaxies with types T — 3 - 5 (63% for the optically complete 
sample and 51% for the CO redshift-restricted sample). 



3. CO observations and analysis 

We carried out CO(1-0) observations with the FCRAO and with 
the IRAM 30m telescope. We observed galaxies with isophotal 
diameters £>25< 100" at the 30m telescope and galaxies with 
£>25> 100" at the 14m FCRAO telescope. In this way we tried to 
optimize the agreement between beam size and optical diameters 
and minimize the fraction of missing flux in the observations 
done with a single pointing. We observed 100 galaxies at the 
FCRAO radio-telescope and 101 at the IRAM 30m telescope. 
In order to check the consistency of the results we observed 15 
galaxies at both telescopes. 

3.1. Observations 

3.1.1. IRAM 30m telescope 

We observed the 12 CO(1-0) line at 1 15 GHz with the IRAM 30- 
meter telescope on Pico Veleta using the dual polarization re- 
ceivers A 100 and B100, together with the 512 X 1 MHz filter- 
banks. The observations were done in wobbler switching mode 
with a wobbler throw of 120" in azimuthal direction. Pointing 
was monitored on nearby quasars every 60 - 90 minutes. The 
integration time on source was typically 0.5 to 1.5 hours. The 
mean system temperatures was 320 K on the T* K scale. All CO 
spectra and intensities are presented on the main beam temper- 
ature scale (r m b) which is defined as r m b = (F e g/B e ff) x 7*£. 
The IRAM forward efficiency, F e s, was 0.95 and the beam ef- 
ficiency, fi e ff, 0.75. The peak intensities of our sources ranged 
between about 10 and 80 mK (r m b). 

Most galaxies were observed at the central position with a 
single pointing. The galaxies with the strongest emission at the 
center (19 galaxies) were mapped along the major axis with a 
spacing of 15", until a root mean square (rms) noise of about 3 
mK was reached for a velocity resolution of 10.6 km s _1 . 

3.1.2. FCRAO 14m telescope 

The observations at the FCRAO were done with the receiver 
SEQUOIA, a 4 x 4 pixel array operating from 85 to 115 GHz. 
We used the so called "beam switching" mode, in which the tele- 
scope switches position between the source and a reference po- 
sition 89.2" apart in azimuth. Two of the pixels in the array al- 
ternated between the ON- and the OFF-position, doubling in this 
way the effective integration time. 

The typical observing time per object was about 2 hours. 
The pointing was checked between the observations of different 
sources using a nearby quasar. The mean system temperatures 
was 380 K (on the T* A scale). All CO spectra and luminosities are 
presented on the main beam temperature scale (r m b) which is de- 
fined as r m b = T^/B e g. The main beam efficiency is B e ff = 0.45. 
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Fig. 1. Basic properties (normalized distributions) of the 
redshift-restricted CO sample (n = 173, full black line), total CO 
sample, excluding possibly interacting objects (« = 234, dashed 
blue line) as well as the optica lly complete sample (n - 712, 
IVerdes-Montenegro et al.l 120051 dotted red line), a) Recession 
velocity, b) optical luminosity, log(L B ) and c) morphological 
type (T(RC3), as in the RC3 catalogue). 



We observed each galaxy with one pointing at their central po- 
sition. The peak intensities of our sources ranged between about 
10 and 80 mK (T mb ). 



Table 2. List of papers from which the CO data have been com- 
piled. 



Bibl. code Article 


Telescope 1 


N 


1 This work 


1,2 


189 


2 Young et al. (1995) 


2, 5 


25 


3 Braine & Combes (1993' 


1 


11 


4 Sautv et al. (2003) 


1,3,5 


on 


5 Elfhae et al. (1996) 


3,4 


15 


6 Sage (1993) 


5 


7 


'The codes for the telescopes are listed in Table[3] 




Table 3. Main parameters of the used radio-telescopes. 




Code Radio-Telescope Diameter ( 


Ihpbw (115 GHz) 


Jy/K 1 


1 IRAM 30m 


21" 


5 


2 FCRAO 14m 


45" 


20 


3 SEST 15m 


43" 


19 


4 Onsala 20m 


33" 


12 


5 NRAO Kitt-Peak 12m 


55" 


32 



The conversion factor from T m b in K to flux in Jy. 



3.1.3. Literature data 

We furthermore searched the literature and found CO(1-0) data 
for 131 objects, 87 of them had not been observed by us. We list 
the references for these data in Table [2] Some galaxies were ob- 
served at several telescopes. Table[3]provides information about 
the telescopes used in the different surveys: the antenna size (col- 
umn 3), the half power beam width (HPBW) (column 4) and the 
conversion factor Jy/K (on the T m b scale) at 115 GHz (column 
5). 



3.2. Data reduction 

The data from both telescopes were reduced in the standard way 
using the CLASS software in the GILDAS packagsfl The data 
reduction consisted in dismissing poor scans, flagging bad chan- 
nels, subtracting a baseline and averaging the spectra for the 
same object and position. In most cases a constant baseline was 
subtracted and only in a few cases the subtraction of a linear 
baseline was required. 



3.3. Spectra and integrated intensities 

The CO(1-0) profiles of the detections and tentative detections 
observed by us at the IRAM 30m and FCRAO 14m are shown in 
Appendix A (Fig. [AT] and Fig. IA.2I respectively). The spectral 
shapes observed are very diverse. Both single and double peaked 
lines are present and the line widths span a wide range. 

The velocity integrated intensity, Iqo = J T mb dv (in K km 
s ), was calculated from the spectra with a velocity resolution of 
10.6 km s _1 for the IRAM spectra, respectively 13.1 km s _1 for 
the FCRAO spectra by summing up all channels with significant 
emission. Its error was calculated as: 



error (7 C o) = cr( W C0 6V C o) 



1/2 



[Kkms- 1 ], 



where <x is the rms noise of the spectrum, Wco is the CO line 
width, and 6Vco is the spectral resolution. For undetected galax- 



http://www.iram.fr/IRAMFR/GILDAS 
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ies, we calculate a 3cr upper limit, assuming a line width of 300 

km s , as: 

/„ , (CO) = 3 x cr(300(5V C o) 1/2 [K km s -1 ]. 
In Table |4] we list the following items: 

1. Entry number in the Catalogue of Isolated Galaxies (CIG). 
An asterisk added to the number means that the detection is 
marginal. In our statistical analysis we treat marginal detec- 
tions as upper limits. 

2. Off. a: RA offset from the center in arcsec. 

3. Off. 6: Declination offset from the center in arcsec. 

4. rms: root mean square noise in mK for a velocity resolution 
of 10.6 kms- 1 (IRAM), respectively 13.1 kms 4 (FCRAO). 

5. Ico'. velocity integrated CO line temperature J T m b dv, in K 
km s _1 , and its error. 

6. Vco: mean velocity of the CO line, in km s _1 . 

7. Wco : zero line- width of the CO spectrum, in km s" 1 . 

8. Tel: radio-telescope code, as listed in Table[3] 

3.4. Comparison between IRAM and FCRAO data 

In order to check the relative calibration between the IRAM 30m 
and the FCRAO 14m telescope and to guarantee that these two 
data sets are comparable, we observed 15 galaxies with both tele- 
scopes. We expect a ratio of the velocity integrated intensities 
of /co-iram/^co-fcrao = 1 for emission homogeneously filling 
the beams, and /co-iram/^co-fcrao = (®fcrao/©iram) 2 = 4.5, 
where ©fcrao and ©iram are the FWHM of the respective 
beams, for a point-like emission. 

Four galaxies (CIG 66, 181, 281 and 330) were detected 
at both telescopes. The ratios of /co-iram/^co-fcrao range be- 
tween 1 . 1 and 2.3, consistent with the value expected for slightly 
concentrated emission. Six galaxies were detected at IRAM, 
but only tentatively detected (CIG 176, CIG 355) or undetected 
(CIG 217, CIG 561, CIG 609, CIG 622) at the FCRAO. The 
lower limit for /co-iram/^co-fcrao in five cases was between 
0.94 and 3.6, consistent with the expected range of values. For 
CIG 217 this value is higher (/co-iram//co-fcrao = 6.1) than 
the theoretical upper limit. Since the detection at IRAM has a 
high signal-to-noise ratio, the most likely reason is an under- 
estimate of the upper limit of the FCRAO data. There is one 
object with a detection at the FCRAO and only a tentative de- 
tection at IRAM (CIG 433), and one with a nondetection at 
IRAM (CIG 268). The ratio of the intensities in both cases is 
^co-iram/^co-fcrao = 0.3, indicating an underestimate in the 
IRAM data. In the remaining three cases, both observations were 
either no detection or tentative detections. 

We conclude that there is very good agreement between the 
detected values at both telescopes, and in most cases (with the 
exception of three galaxies) also for objects only detected at one 
telescope. This gives us confidence that the calibration of the two 
data sets is consistent. 



3.5. Calculation of the molecular gas mass 

The molecular gas mass (Mh,) is calculated using a Galactic 
con version factor of NfH 2 )/Ico = 2.0 x 10 2() cm ~ 2 (K km s~V 
(e.g. lDickman et al.1 ll986) yielding: 



the observations in arcsec 2 (i.e. Q = 1.1 30^ for a single point- 
ing with a Gaussian beam of FWHM 0#). We do not include the 
mass of heavy metal (mostly helium) in the molecular gas mass. 

Most of our objects were observed at the central position in 
a single pointing since the mapping of the entire galaxy would 
have been too time-consuming. We therefore might have missed 
part of the CO emission for galaxies where the emission is more 
extended than the beam. This fraction depends on the galaxy 
size, inclination and on the telescope beams. It is thus necessary 
to correct for this loss, and we do this by extrapolating Mh, ob- 
served in the central beam to the total mass in the galaxy. In the 
next subsection we explain how we carried out this correction. 

3.5.1. Aperture correction 

In order to apply an aperture correction, we need to predict the 
distribution of the CO emis sion. CO maps of nearby spiral galax- 
ies dNishivam a et al.ll200UlRegan et alfeooil ILerov et al.l l2008) 



have shown that the radial distribution of Ico(r) in galaxies can 
be well described by an exponential function with a scale length 



Ico(r) = lo exp(-r/r e ). 



(2) 



M H ,[M S ] = 75/ co D 2 Q, 



(1) 



where Iqo is the velocity integrated CO line intensity in K 
km s _1 , D is the distance in Mpc and D. is the area covered by 



The CO scale length, r e , is w ell correlated and similar to the 
optica l exponential scale length (Reg an et alJl200lt iLerov et all 
2008). It also correlates, although l ess tightly, with the optical 
radius at the 25mag isophote, r25. ILerov et ail d2008l) derived 
for spiral galaxies from the THINGS survey a mean value of 
oe = r e /r25 = 0.2. We derived the same mean value for a from 
the data of lRegan et al.l (|2001|) for 15 spiral galaxies observed in 
the BIMA Survey of Nearby Galaxies (BIMA-SONG) and from 
the data of Nishiyama & Nakal d2001f) for 25 spiral galaxies ob- 
se rved with the Nobeyama 45m telescope. We also used the data 
of I Young et al.l d 19951) . who studied the molecular gas content 
and distribution in a sample of 300 nearby galaxies, to derive 
the r e . They found a mean ratio between the effective CO diam- 
eter, Deo, the diameter within which 70% of the CO emission 
is situated, and the optical diameter of DC0/D25 = 0.5. For an 
exponential distribution one can derive that Deo x 0.5 = 2.5r e . 
Thus, their data also yield r e /r 25 = Dc £° x lf 5 = 53 = 0.2. 

Finally, we use the data of the 19 galaxies (all of them with 
morphological type T > 2) mapped along the major axis with the 
IRAM 30m telescope as a further test. Although our data is not 
sufficiently detailed to fit the radial distribution (we have only 
3-5 detected spectra along the major axis), we can use it to test 
whether (i) an exponential distribution is a reasonable descrip- 
tion of the CO distribution, (ii) the scale length derived by other 
studies is in agreement with our data, and (iii) the predictions for 
the extrapolated Mh, are in agreement with our mapped values. 

For the first two tests, we fitted an exponential distribution 
independently to each side of the CO distribution along the ma- 
jor axis. From the 38 resulting fits, there were only six cases 
where no exponential fit could be applied within the error bars. 
In five cases we derived a > 0.4, in seven cases 0.4 < a < 0.25, 
in three cases a < 0.15 and in 17 cases, the majority, 0.15 < a < 
0.25. Thus, our data are in general consistent with the value of a 
found by other studies. 

In summary, we conclude that an exponential distribution of 
the molecular gas distribution with a = r e /r25 = 0.2 is a good 
approximation based on the CO maps for nearby spiral galaxies 
available up to date and also consistent with our data. We adopt 
the same value for both spiral galaxies and early type galax- 
ies. In early-type galaxies, the molecular gas extent is much less 



6 U. Lisenfeld et al.: The AMIGA sample of isolated galaxies. X 

Table 4. Velocity integrated 12 CO(1-0) line intensities, mean velocity and line widths for the CIG galaxies observed by us. 



CIG 


Off. a 


Off. S 


rms 


Ico 


Vco 


W co 


Tel 




["] 


["] 


[mK] 


[Kkms- 1 ] 


[kms- 1 ] 


[krns" 1 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


10* 








3.47 


0.50 ±0.15 


4980 


180 


1 


27 








4.21 


3.54 ± 0.23 


4586 


287 


1 


27 


13 


8 


7.21 


< 1.22 






1 


29* 








5.24 


0.49 ±0.16 


4162 


92 


1 



The full table is available in electronic form at the CDS and from http://amiga.iaa.es 



known. However, for our study, this uncertainty is not important 
because the number of objects with T < is low (n = 23, with 
eight detections) and we focus on our results and conclusions on 
spiral galaxies, in particular of types 7 = 3 — 5, which dominate 
the sample. 

We now use these results to calculate the aperture correc- 
tion which we define as the ratio between the total (extrapolated) 
molecular gas mass Mh,, and the molecular gas mass in the cen- 
tral pointing, M H ,,center, 



Zap = Mh,/M H ,,< 



(3) 



The total molecular gas mass is calculated by spatially inte- 
grating Ico{f) from Eq.|2]and using Eq.[TJ This yields : 



M H ,[M ] = 75D 2 I{r)2ndr 



Jo 

T 

Jo 



(4) 



75D 2 I /oexp(-r/r e )27T(/r = 75D 2 / 27rr 2 . 



Similarly, we calculate Me 2l center by convolving the exponen- 
tial CO intensity distribution with a Gaussian beam. This yields: 



Mi 



H2, center L 



r [M Q ] =75D 2 4/ g dx f Q °° dy 
e X p(-ln(2)[(^)^(^) 2 ])e XP (-4^)^ 



(5) 



where i is the inclination of the disk. The integration of Eq.[5]is 
carried out numerically. 

Thus, the correction factor, / ap =MH,/MH 2;C enter, depends on 
the ratio of the scale length and the beam size, r e /&B, as well as 
the galaxy inclination i. Fig. [2] shows the distribution of / ap for 
the galaxies in our sample. The correction factors are generally 
low: 81% of the galaxies have / ap < 2, and 92% / ap < 3. Only 
nine galaxies have a correction factor above 5. All of them are 
nearby (v < 1000 km s~') galaxies with a large angular size 
(between 4' and 20') that are not included in our redshift-limited 
sample. 

In order to carry out test (iii) we compared for the 19 galax- 
ies, that were mapped along their major axis with the IRAM 
30m telescope, the extrapolated Mh, to the mapped molecular 
gas mass which we extrapolated to the mass in the entire disk by 
assuming azimuthal symmetry (MH,, m ap,extra)- Fig.[3]presents the 
ratio of these two masses as a function of D25 and 2. For most ob- 
jects, both masses agree reasonably well, with ratios ranging be- 
tween 0.5 and 1.3. There are only two outliers, with ratios around 
2. These galaxies (CIG 84 and CIG 28) have a very flat CO dis- 
tribution along the major axis so that our extrapolation, assum- 
ing an exponentially decreasing distribution, underestimates the 
true amount of molecular gas. No trend with neither D25 nor the 



2 3 

aperture correction factor f 



Fig. 2. Histogram of the aperture correction factor, / ap , for galax- 
ies observed by us with the IRAM 30m telescope (upper panel), 
with FCRAO (middle panel) and for galaxies taken from the lit- 
erature (lower panel). 



inclination is seen, showing that no apparent bias is introduced 
by the aperture correction. The mean mass ratio is 1.0, with a 
standard deviation of 0.3. 



4. Results 

4.1. Molecular gas content 

In Table [5] we list Mh, for the individual galaxies. If observa- 
tions from different references were available for a single galaxy, 
we inspected the spectra and discarded those of poorer quality. 
In case of similarly good data we gave preference to mapped 
observations or to observations from the telescope with a larger 
beam in order to avoid flux loss. For the 44 objects that had data 
from the literature and were also observed by us, our data was 
in general of better quality, with the exception of three objects 
(CIG 512, CIG 604, CIG 626). The columns in Table |5]are: 

1. CIG: Entry number in the Catalogue of Isolated Galaxies 
(CIG). 

2. log(Mn 2 ) center: Decimal logarithm of Mh, towards the cen- 
tral pointing in solar masses, derived from Eq. [TJ An asterix 
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Fig. 3. Ratio between the molecular gas mass, mapped along the 
major axis with the IRAM 30m telescope, and extrapolated to 
the entire disk assuming azimuthal symmetry (MH 2>m ap,extra) and 
the molecular gas mass extrapolated from the central pointing as 
explained in Sect. 3.5.1 (Mh 2 ) as a function of diameter (a) and 
galaxy inclination (b). 

Table 5. Molecular gas mass 



CIG 




log(M H2 )[M ] 


Tel. 


Ref. 




center 


mapped 


extrapol. 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


9.47 





9.58 


3 


4 


4 


8.96 


9.08 


9.26 


2 


2 


6 


<8.14 





<8.14 


5 


4 


10 


7.88* 





8.13* 


1 


1 



The full table is available in electronic form at the CDS and from 
http://amiga.iaa.es 



denotes tentative detections that are treated as upper limits in 
the statistical analysis. 

3. log(Mn 2 ) mapped: Decimal logarithm of the mapped Mh 2 in 
solar masses, calculated in the following way: for the data 
from the literature, the angular separation between the indi- 
vidual pointings was always larger than the beam size so that 
the total Mh 2 could be calculated as the sum of the individ- 
ual pointings. For our own observations with the IRAM 30m 
telescope, the spacing between the individual pointing was 
15", which is smaller than the FWHM of the beam (21") so 
that in this case we had to take the overlap of the individual 
pointings into account. We calculated the mapped Mh, from 
Eq. [TJ where Iqo is taken as the mean value of the different 
pointings, and Q is the total area covered by the mapping, 
approximated as 21" x 36", 21" x51", 21" x66", 21" X 81" 
for 2, 3, 4, and 5 pointings, respectively. 

4. log(Mn 2 ) extrapol.: Total (extrapolated) Mh, in solar masses, 
calculated as described in Section 3.5.1. An asterisk denotes 
tentative detections that are treated as upper limits in the sta- 
tistical analysis. 

5. Tel.: Radio-telescope code, as in Table[3] 

6. Ref.: Bibliographic code, as in Table [2] 

Fig. [4] shows the distribution of log(MH,fl for detections in 
the total and redshift-limited CO samples and for nondetections 
in the latter. The distributions in the total and in the redshift- 
limited samples are very similar. The mean Mh, of the redshift- 
limited sample is log(Mn 2 ) = 8.30 ± 0.08 (see Table [§}, calcu- 



2 Here, and in the following, we always use to the extrapolated 
molecular gas mass in our analysis and denote it Mh, for simplicity. 



20 



i 



log(M H2 ) [MJ 



10 



Fig. 4. Histogram of the molecular gas mass, Mh 2 . Galaxies 
flagged as potentially interacting were excluded. The full line 
corresponds to all detections (n - 131), the gray filled histogram 
shows detections in the redshift limited sample (n = 89) and 
the dashed lines shows tentative or nondetections in the redshift- 
limited sample (n = 84). 
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Morphological type T(RC3) 

Fig. 5. The molecular gas mass for the redshift-limited sam- 
ple as a function of morphological type. Triangles denote the 
mean value and its error for a range of morphological types and 
squares the median value, as listed in Table [6] The error bars 
in the x-direction denote the range of morphological types over 
which the mean and median have been taken. 



lated with the package ASURV that takes into account the upper 
limits. Here and throughout the paper we use the AS URV0 pack- 
age which applies survival analysis in the presence of upper or 
lower limits and calculates the mean value based on the Kaplan- 
Meier estimator. 



3 Astronomy Survival Analysis (ASURV) Rev. 1.1 is a general- 
ized statistical package th at i mplements t he m ethods presented by 
Feigelson & Nelson ( 1985) and Isobe et al. ( 1986), and is described in 
detail in llsobe et al] fi99Ch and lLavallev et al.Ul992h . 
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Fig. [5] shows the distribution of the molecular gas mass as a 
function of the morphological type. The mean and median values 
are listed in Table[6] Our sample is dominated by spiral galaxies 
of type T — 3 -5 (Sb-Sc). Not only is the total number of objects 
greatest in this range, but also the detection rate. Therefore, we 
can derive the most reliable results for these types. Both for ear- 
lier and for later types, the detection rates are very low, making 
a detailed analysis difficult. The molecular gas mass is largest 
for spiral galaxies of T = 3-5, and decreases both for earlier 
and later types. There are eight early-type galaxies, of type SO 
and SOa, with detections in Mh, , and five of them have unusu- 
ally high molecular gas masses in the range of those for spiral 
galaxies (CIG 332, CIG 481, CIG 498, CIG 733 and CIG 1015). 

4.2. Relation of Mh 2 to other parameters 

In the following we investigate the relations between Mh 2 and 
Lb, D25, Lpm, Lk and Mhi. The first two quantities (Lb, D25) 
were chosen because they are in general available for any galaxy 
and are therefore useful to predict the expected Mh 2 . Lfir is very 
closely related to Mh, because of their common relation to SF. 
Lk is dominated by the emission of low-mass stars which are 
the result of the long-term SF history of an object and deter- 
mine the gravitational potential which influences the SF activity. 
Finally, we compare Mh, to Mhi in order to derive the molecular- 
to-atomic gas mass ratio as a function of morphological type. 

The results of the regression analysis as well as the Spearman 
rho correlation coefficient derived in this section are listed in 
Table [7] and the mean and median values of the ratios for dif- 
ferent subsamples are in Table [6] We use the package ASURV 
to calculate the bisector regre ssion line applying the Schmitt's 
binning method (Schmittl ll985l) as the only method offered by 
ASURV able to deal with censored data in both the dependent 
and the independent variably 

4.2.1. Optical luminosity 

Fig. [6] shows the relation between Mh, and Lb, for the the en- 
tire redshift-limited sample and for subsamples of different mor- 
phological types, together with the the best-fit bisector regres- 
sion lines. The high number of upper limits for galaxies of type 
T = (-5) -2 and T =6-10 impedes the calculation of a reliable 
regression line for these individual subgroups. However, we note 
that there is no apparent deviation from the mean best-fit regres- 
sion for all morphological types for these two groups. A slope 
considerably larger than 1 is fo und for the corre lation between 
Mh, and Lb, in agreement with lPerea et al.l d 19971) . This means 
that the ratio Mh,/Lb tends to increase with Lb. Even though 
the scatter around the best-fit is large, so that no strong corre- 
lation exists between Mh,/Lb and Lb, we confirm the variation 
of Mh,/Lb by finding that the mean values is indeed higher for 
high Lb than for low Lb: for galaxies (all morphological types) 
with Lb< 10 10 L Q (n = 123 galaxies, n = 68 upper limits) we 
obtain M H ,/L B =-1.36 ± 0.05 and for galaxies with L B > 10 10 L 
(n/n =50/11) Mh 2 /L b = -1.10 ± 0.06. The corresponding num- 
bers for the T = 3 - 5 subsample are -1.16 + 0.06 (L B < 10 10 
L , n/n =48/16) and -1.04 ± 0.05 (L B > 10 10 L , n/n =40/5). 
This trend has to be taken into account when using this ratio as 
an indicator of an enhancement of Mh, . 

The ratio between Mh, and Lb is shown in Fig. [7] as a func- 
tion of morphological type. The values are listed in Table|6] The 

4 The Schmitt method was partially reimplemented and wrapped into 
Python. It can be found at http://amiga.iaa.es/software/python-asurv 
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Fig. 6. M H , vs. L B for the redshift-limited CO sample (n = 173), 
including all morphological types, and for groups of different 
morphological types. The dashed line gives the best fit bisector 
(derived with ASURV) for all morphological types and the solid 
line is the best fit for the T = 3 - 5 sample. 
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Morphological type T(RC3) 

Fig. 7. The ratio between Lb and Mh, as a function of morpho- 
logical type. Triangles denote the mean value and its error for 
a range of morphological types and squares denote the median 
values, as listed in Table [6] The error bars in the x-direction de- 
note the range of morphological types over which the mean and 
median have been taken. 



ratio is highest for galaxies of type T - 3 - 5. It is however re- 
markable that galaxies of type SO and SOa with detections in CO 
have much higher values of Mh,/Lb, in the range of spiral galax- 
ies of type T = 3 - 5, than early-type galaxies with nondetection 
in CO. 

In order to check whether the extrapolation of the molec- 
ular gas mass to the entire disk has introduced any biases, we 
show in Fig. |8]the relation between Mh, and Lb only for galax- 
ies with / a p < 1.5. No significant difference compared to the 
entire redshift-limited CO sample can be seen, and the best-fit 
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Fig. 8. Mh, vs. Lb for galaxies in the redshift-limited CO sample 
with an aperture correction factor of less than / ap < 1.5 (n — 78). 
The solid line is the best fit bisector line derived with ASURV 
for this restricted sample, and the dashed lines gives, for com- 
parison, the fit for the entire redshift-limited CO sample. 



regression coefficients are the same within the errors (Table [7]). 
We found this good agreement between entire redshift-limited 
CO sample and the subsample of galaxies with / ap < 1 .5 for all 
correlations studied. Therefore, in the following subsections we 
do not show the f ap < 1 .5 correlation separately, but we list in 
Tables [6] and Q the mean values and the regression coefficients 
for this subsample. 
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Fig. 9. Mh, vs. square of the optical isophotal diameter, 
, for the redshift-limited CO sample (« = 173), including all 
morphological types, and for groups of different morphological 
types. The dashed line gives the best fit bisector (derived with 
ASURV) for all morphological types and the solid line is the 
best fit for the T = 3 - 5 sample (see Table|7]). 



4.2.2. Optical isophotal diameter 

Fig. [9] shows the relation between Mh, and for the entire 
redshift-limited sample and for subsamples of different morpho- 
logical types, together with the best-fit bisector regression lines. 
For early-type galaxies (T < 0) only a very poor correlation is 
visible, whereas galaxies of type T = 6-10 seem to follow the 
same correlation as those of type T = 3 — 5. 

The correlation between Mh, and D^ 5 , has the lowest corre- 
lation coefficient (r ~ 0.5) among those considered in this paper 
(Table 17). Although the bisector slope is formally larger than 1, 
we do not find a variation of MuJD^ with increasing optical di- 
ameter. This shows that for this poor correlation the regression 
slope has to be taken with caution. 

Fig. [10]shows the ratio MnJEfc* as a function of morpholog- 
ical type. The values are listed in Table|6] The ratio is highest for 
spiral galaxies of type T = 3 — 5. Similarly to Mh,/Lb, SO and 
SOa galaxies with detections in CO have high values of MuJD^, 



similar to those of spiral galaxies of type T 



Mh,/£>L of the nondetections is much lower. 



3-5, whereas 



4.2.3. Luminosity in the K-band 

Fig. [TT] shows the relation between Mh, and Lk, for the entire 
redshift-limited sample and for subsamples of different morpho- 
logical types, together with the the best-fit bisector regression 
lines. The relation is close to linear. The correlation is very good 
(r = 0.73) for spiral galaxies of type T = 3 — 5. The distribution 
of the emission from spiral galaxies of later types (7 = 6 — 10) 
is consistent with this correlation. However, early-type galaxies 
(T < 0) show a very poor correlation. Only the objects with CO 
detections follow the same correlation as spiral galaxies whereas 
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Fig. 10. The ratio between Mh 2 and D^ 5 for the redshift-limted 
sample as a function of morphological type. Triangles denote the 
mean value and its error for a range of morphological types and 
squares the median values, as listed in Table [6] The error bars 
in the x-direction show the range of morphological types over 
which the mean and median have been taken. 



most galaxies with CO nondetections have upper limits for Mh, 
that lie considerably below it. 

The ratio Mh 2 /Lk is shown in Fig. [12] It is lowest for early- 
type galaxies (up to T = 2), and increases for later types by a fac- 
tor 3-5. From T = 3 on, the ratio is approximately constant and 
does not show the decrease it shows for Mh,/Lb and Mh 2 /D?,. 
Early-type (E+S0s) galaxies have a lower mean molecular gas 
mass per Lk than spiral galaxies. This is caused by galaxies not 
detected in CO, whereas, as seen before, SO and SOa detected in 
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Table 6. Mean and median values. 



Sample 


< loa(Mu,)> 


< log(^S) > 


< log(^P) > 

25 , 


< log(^) > 


< log(^S) > 


< loe(^ffi-) > 

& v M H2 ' 


<log(SFE)> 




[M Q ] 


[Mo/L ] 


[M /kpc 2 ] 


[Mo/L^o] 




[L e /M ] 


[yr '] 




median 


median 


median 


median 


median 


median 


median 




n/n up 


w/«up 


«/«up 


n/n up 


n/n up 


n/n up 


«/n up 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


E-Im (r=(-5) - 10) 


8.30±0.08 


-1.28+0.04 


6.00+0.05 


-1.87+0.05 


-0.94+ 0.05 


0.72+0.03 


-8.94 




8.61 


-1-10 


6.22 


-1.58 


-0.73 


0.64 


-9.02 




173/79 


173/79 


173/79 


157/65 


153/63 


97/22 


97/22 


E-Im(/ ap < 1.5) 


8.30+0.12 


-1.29+0.06 


6.06+0.07 


-1.80+0.06 


-0.96+0.08 


0.80+0.06 


-8.86 


8.51 


-1.07 


6.29 


-1.55 


-0.76 


0.69 


-8.97 




76/42 


76/42 


76/42 


65/33 


66/33 


39/12 


39/12 


E-SOa (T=(-5) - 0) 


8.03±0.12 


-1.45+0.11 


5.79+0.13 


-2.32+0.10 


-0.46+0.21 


0.95+0.08 


-8.71 




8.35 


-1.16 


6.27 


-2.06 


-0.25 


0.94 


-8.72 




23/15 


23/15 


23/15 


22/15 


8/3 


6/0 


6/0 


Sa- Sab (7=1-2) 


8.49+0.18 


-1.41+0.19 


5.76+0.20 


-2.27+0.19 


-0.58+0.21 


0.69+0.06 


-8.97 




O.OO 


-uz 


j.oy 


-Z.Uo 


-U.53 


U.OJ 


-9.03 




9/5 


9/5 


9/5 


9/5 


8/4 


6/3 


6/3 


Sb-Sc (T= 3-5) 


8.79+0.06 


-1.10+0.04 


6.26+0.05 


-1.66+0.04 


-0.72 ±0.06 


0.63+0.03 


-9.03 




8.87 


-1.01 


6.37 


-1.58 


-0.60 


0.62 


-9.04 




88/21 


88/21 


88/21 


87/20 


85/19 


68/10 


68/10 


Scd-Im(T= 6-10) 


7.78+0.15 


-1.41+0.05 


5.79+0.05 


-1.62+0.06 


-1.38+0.07 


0.99+0.09 


-8.67 




8.26 


-1.15 


6.01 


-1.37 


-1.07 


0.63 


-9.03 




53/38 


53/38 


53/38 


39/25 


52/37 


17/9 


17/9 



Mean value and its error calculated with the program ASURV for different subsamples, and (below) median value (calculated treating the upper 
limits as detections). For the mean and median value of log(M H ,/LFiR) and log(M H2 /M m ) only galaxies with detection in L FIR or M m , respectively, 
were taken into account, since ASURV cannot handle both upper and lower limits. Below: Total number of galaxies and number of upper limits in 
Mh 2 taken into account for the means and median. 
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Fig. 11. The relation between Mh, and Lk for the redshift- 
limited CO sample (n = 173), including all morphological types, 
and for groups of different morphological types. The dashed line 
gives the best fit bisector (derived with ASURV) for all morpho- 
logical types and the solid line is the best fit for the T — 3 - 5 
sample (see Table|7]). 



CO have higher values of MuJL^, in the same range as spiral 
galaxies. 
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Fig. 12. The ratio between Mh, and Lk as a function of morpho- 
logical type. Triangles denote the mean value and its error for 
a range of morphological types and squares denote the median 
values, calculated treating upper limits as detections, as listed 
in Table [6] The error bars in the x-direction denote the range of 
morphological types over which the mean and median have been 
taken. 



4.2.4. FIR luminosity, SF rate and efficiency 

A good correlation is kn own to exist between Mh, and Lfir (e.g. 
lYoung & Scovillell 19911) because both quantities are directly re- 
lated to SF: the molecular gas as the fuel for SF and Lfir as a 
tracer for SF based on the heating of the dust by newly born stars. 
Fig. [13] shows the relation between Mh, and Lfir, for the entire 
redshift-limited sample and for subsamples of different morpho- 
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Table 7. Correlation analysis 



Sample 


n 


n up 


slope 


intercept 


slope 


intercept 


corr. coeff 








(bisector) 


(bisector) 


(M H , dep.) 


(M H , dep.) 


r 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Mr 2 vs. L B 


E-Im (T=(-5) - 10) 


173 


79 


1.45± 0.08 


-5.61±0.77 


1.12+ 0.08 


-2.43+0.83 


0.66 


E-Im(/ ap < 1.5)* 


76 


42 


1.32±0.10 


-4.27±0.94 


1.06+0.10 


-1.79+0.98 


0.63 


Sb-Sc (T = 3-5) 


88 


21 


1.41±0.10 


-5.12±0.97 


1.06+0.10 


-1.69+1.02 


0.65 


M H , vs. Di s 


E-Im (L=(-5) - 10) 


173 


79 


1.41±0.05 


5.10±0.50 


0.88+0.09 


6.36+0.23 


0.52 


E-Im(/ ap < 1.5)* 


76 


42 


1.32±0.08 


5.47±0.71 


0.91+0.11 


6.38+0.26 


0.47 


Sb-Sc (3-5) 


88 


21 


1.31±0.09 


5.50±0.93 


0.79+0.10 


6.82+0.27 


0.53 


M H , vs. L K 


E-Im (L=(-5) - 10) 


157 


65 


1.05±0.07 


-2.27+0.72 


0.72+0.07 


1.13+0.74 


0.64 


E-Im (/ ap < 1.5)* 


65 


33 


0.99+0.09 


-1.54+0.96 


0.69+0.09 


1.41+0.95 


0.66 


Sb-Sc (T = 3 - 5) 


87 


20 


1.18±0.07 


-3.49+0.78 


0.95+0.08 


-1.14+0.85 


0.74 


M H , vs. L FIR 


E-Im (L=(-5) - 10) 


172 


97 


1.16±0.08 


-2.14+0.72 


0.98+0.06 


-0.46+0.61 


0.80 


E-Im(/ ap < 1.5)* 


75 


48 


1.15±0.13 


-2.09+1.23 


0.89+0.11 


0.28+1.04 


0.71 


Sb-Sc (T = 3 - 5) 


88 


30 


1.04±0.06 


-1.00+0.55 


0.90+0.06 


0.29+ 0.58 


0.82 


Zsfr vs. 


E-Im (L=(-5) - 10) 


172 


97 


0.89±0.07 


-2.78+0.20 


0.82+0.07 


2.25+0.20 


0.65 


E-Im(/ ap < 1.5)* 


75 


48 


1.00±0.14 


-2.77+0.38 


0.57+0.10 


1.53+0.31 


0.56 


Sb-Sc {T = 3 - 5) 


88 


30 


0.88±0.07 


-2.82+0.21 


0.84+0.08 


2.39+ 0.22 


0.71 



The entries are: Column 1: Subsample considered. Column 2: Total number of galaxies in the respective samples. Column 3: Number of galaxies 
with upper limits in Mh, and/or Lfir. Column 4: Bisector slope and its error of the best-fit regression line derived with the Schmitt binning method 
in the ASURV package. The slope and intercept are defined as log(A) = intercept + log(B) x slope, where A is M H , or Z H2 and B is L B , L K , 
Lfir, (D25) 2 or Ssfr, respectively. Column 5: Bisector intercept and its error. Column 6 and 7: Slope and intercept and their errors of the best-fit 
regression line derived with the Schmitt binning method in the ASURV package adopting L B , L F ir , L K D^ 5 or E SF r as independent variable. 
Column 8: Spearman's rho correlation coefficient, calculated with ASURV. 
* Only galaxies for which the aperture correction factor, / ap (see Sect. 3.5.1), is less than 1.5. 



For early (T < 0) and late (T — 6 - 10) type galaxies, the dis- 
tribution of Mh, vs. Lfir is consistent with the correlation found 
for galaxies of type T = 3 - 5. 

Fig. [14] shows the ratio Lfir/Mh, as a function of morpho- 
logical type. Early-type galaxies have a higher value, possibly 
due to a large fraction of their FIR emission not being heated 
by young stars but by the general interstellar radiation field. For 
late-type spirals (T > 6) the mean ratio increases again, but the 
low number of values and detections make any firm conclusion 
difficult. 

These val ues and trends for Lf ir /Mh, are consistent with ear- 
lier results of I Young et alj (Q~996) who studied them in a sam- 
ple of 120 galaxies included in the FCRAO survey. They used a 
higher value of N(H 2 )/I C o = 2.8 x 10 20 cirT 2 (Kkms^'r 1 , but 
also define the FIR luminosity in the range between 1 - 500 fim 
which is, according to their prescription, between ~ 0.1 — 0.2 
dex higher. Since both differences roughly compensate, their val- 
ues of Lir/Mh 2 are comparable to our values Lfir/Mh,. They 
find very similar values of Lfir/Mh, for different morpholog- 
ical types as we do. They obtain value between 0.55 ± 0.08 
and 0.61 + 0.06 for morphological types T — 3 - 5 (n = 45), 
higher values (between 0.70 + 0.13 and 1.53 ± 0.24) for later 
type galaxies (T = 6 - 10, n — 19), and, for earlier spiral types 
(T = 1 - 2, n = 14), mean values (0.65 ± 0.20 and 0.53 + 0.09) 
that are similar to galaxies of type T = 3 — 5. Their sample only 
includes three galaxies of type SO-SOa and no elliptical galaxies, 
so that we cannot compare earlier types. 

Lfir is a good tracer for the star formation rate (SFR) due to 
two reasons: (i) young stars are formed in dense regions where 
dust opacities are high and (ii) the dust extinction curve peaks 
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Fig. 13. Mh, vs. Lfir for for the redshift-limited CO sample 
(n = 173), including all morphological types, and for groups 
of different morphological types. The dashed line gives the best 
fit bisector (derived with ASURV) for all morphological types 
and the solid line is the best fit for the T — 3 - 5 sample (see 
Tabled). 



logical types, together with the best-fit bisector regression lines. 
We find a good correlations with a roughly linear slope (Table|7}. 
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Fig. 14. The ratio between Lfir and Mh 2 as a function of mor- 
phological type. Galaxies with upper limits at both wavelengths 
are excluded. Triangles denote the mean value and its error for 
a range of morphological types and square signs denote the me- 
dian values, as listed in Table[6] The error bars in the x-direction 
show the range of morphological types over which the mean and 
median have been taken. 



in the ultraviolet so that the radiation from young, massive stars 
is preferentially absorbed. Therefore, Lfir as a good SFR tracer 
especially for actively SF galaxies. It has to be used with some 
caution in galaxies with a low SFR where dust heated from old 
stars can contribute to Lfir, or i n galaxies w ith a low metallicity 
and thus a low dust opacity (e.g. Bell 2003). 

Keeping the se limitations in mind, we use the formula of 
iKennicuttl (1 19981) to calculate the SFR: 



SFR(M Q yr ) 



^xlO^Wergyr- 1 ), 



(6) 



where Lir is the total FIR luminosity in the range 8 -1000 fim. 
This formula assumes a Salpet er Initial Mass Function. We con- 
vert this to a val ue based on thelKroupal d2001l) IMF by deviding 
by a factor 1.59 dLerov et al.ll2008l) . In our analys is we use Lfir , 
calculated following the formula given bv lHelouet all (1988). 
which estimates the FIR emission in the wavelength ran ge of 
42-12 2.5 yum. We estimate Lir from Lfir using the result of Bell 
(2003?) that on average Lir ~ 2x Lfir for a heterogenous sample 
of normal and starbursting galaxies. Adopting this factor and the 
conversion to the Kroupa IMF we calculate the SFR from the 
Lfir as: 



SFR(M yr 1 ) = 4.5 x 2 x — lO^LpiR^rgyi-- 1 ) (7) 
= 2.2 x 1O- 1O L FIR (L ). 



We derive an average SFR for galaxies of type T — 3 - 5 
of ~ 0.7 M yr~\ The ratio Lfir/Mh 2 is proportional to the star 
formation efficiency (SFE), defined here as the ratio between 
SFR and Mh 2 . The mean value of the SFE for galaxies of type 
T = 3 - 5 is 10~ 9 yr 1 (Table©. 
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Fig. 15. The relation between Mh, and Mhi for groups of differ- 
ent morphological types. 




5 
Morphological type T(RC3) 

Fig. 16. The Mh,/Mhi ratio as a function of the morphological 
type. Triangles denote the mean value and its error for a range 
of morphological types and squares denote the median values, 
as listed in Table [6] The error bars in the x-direction denote the 
range of morphological types over which the mean and median 
have been taken. 



4.2.5. Atomic gas mass 

Fig. Q3] shows the relation between Mh, and Mhi for differ- 
ent groups of morphological types. The correlation with Mhi is 
much poorer than the other correlations considered here, with 
a Spearman correlation coefficients of r = 0.44, 0.44, 0.29, 0.57 
for the entire redshift-limited sample, T - (-5) — 2, T = 3 - 5 
and T = 6 - 10, respectively. In most galaxies Mhi is higher than 
Mh, • Some differences can be seen for the different morpholog- 
ical types: Early type galaxies (T = — 5 - 2) have a high number 
of upper limits in both Mh, and Mhi- The upper limits mainly 
populate the low Mhi and Mh, part of the diagram, whereas the 
detections have values of Mh, and Mhi comparable to those of 
spiral galaxies. Late-type galaxies (T = 6 - 10) are shifted to- 
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wards higher M m ( > 10 9 M Q ) but lower ( < 10 9 M Q ) M Hl with 
a large number of upper limits in the latter. Galaxies of morpho- 
logical types T = 3 - 5 have both high atomic gas and molecular 
gas masses ( i 10 9 M Q ) with a low number of upper limits in 
Mh, and almost none in Mm. 

Fig. [16] displays the ratio Mh 2 /Mhi as a function of mor- 
phological type, showing a strong variation with T. The high- 
est values are found for early-type galaxies, up to T = (albeit 
with a high uncertainty due to the high number of upper limits). 
For later types, Mh,/Mhi is decreasing strongly. The mean ra- 
tio MuJMhi is significantly lower than 1 for all morphological 
types. 

4.3. Expected molecular gas content in a galaxy 

We showed in Sect. 4.2 that good correlations exist between Mh, 
and other parameters of a galaxy (Lb, Lk and Lfir) and a some- 
what poorer correlation with D^. In Table [7] we list, apart from 
the linear regression parameters, also the Spearman's rho corre- 
lation coefficient and in Table [6] the different ratios depending 
on the morphological types. All these relations can be used to 
study differences in the molecular gas content of other sample, 
like e.g. interacting galaxies, with respect to isolated galaxies. 

The best correlation exists between Mh, and Lfir('" ^ 0.8). 
Lfir is thus a very reliable way for predicting the expected 
molecular gas content in a galaxy. However, Lfir might not be 
a good parameter when searching for variations of Mh, in in- 
teracting galaxies because Lfir itself, tracing the SFR, is easily 
affected in such an environment. A very good, and roughly lin- 
ear, correlation also exists between Mh, and Lk for galaxies of 
type T — 3 - 5 (r = 0.73). The luminosity in the K-band as a 
measure of the total stellar mass is less affected by recent events 
than Lfir or Lb and thus a good normalization parameter when 
searching for changes in Mh, . This correlation is, however, poor 
for early type galaxies (T < 2), where the undetected objects 
have ratios Mh,/Lk well below the values for spiral galaxies with 
T > 3 (see FigQT]and Tabled so that for those types Lk is not 
recommened as a measure of the expected Mh,. There is also 
a good correlation between Mh, and Lb (r ~ 0.65). The ratio 
Mh,/Lb is however not constant, but increases with Lb which 
has to taken into account in any comparison. Finally, the poorest 
correlation (r ~ 0.5) exists with the isophotal diameter which is 
not a very reliable parameter for predicting Mh, . 

Apart from using the ratio of Mh 2 to another parameter, we 
can use the correlations, defined by the linear regression param- 
eters listed in Table [7] for predicting the expected Mh, and de- 
termine whether a deficiency or an excess of Mh, exists in an 
object in comparison to isolated galaxies. 

We define the d eficiency in Mh, in an anal ogous way to the 
definition in Mm of lHavnes & GiovaneliU (e.g. [19841) : 

Def (Mh, ) = log(M H , .predicted) - log(M H2 ,observed)- (8) 

Note that in this definition, a positive value of Def(Mn,) 
means a deficiency and a negative value means an excess. We 
can derive the predicted Mh,, MH 2;Pre dicted, from any of the pa- 
rameters, X, (L B , L>2 5 , L K or Lfir) as: 

log(M H ,, predicted) = intercept + slope X log(X), (9) 

where intercept and slope are the parameters of the best fit listed 
in Table [7] These definitions allow us to directly take the effect 
of a nonlinearity of a correlation into account. 

In Table [7J we have given the values for two types of regres- 
sions: The bisector fit and the fit obtained by taking Mh, as the 



dependent variable, 0(Y|X), and minimizing the distance of the 
Mh, measurements from the best-fit regression line. The regres- 
sion parameters are different for both methods because of the 
scatter in the data and the difference is larger for poorer cor- 
relations. The 0(Y|X) fit is the appropiate regression for pre- 
dicting Mh, from Lb, whereas the bisector regression is the best 
estimate for the un derlying correlation between two parameters 
dlsobe et al.ll 19901) . For a sample covering the same luminosity 
range, 0(Y|X) is the best way for predicting the expected Mh 2 . 
For samples with a different luminosity range, however, the bi- 
sector fit is better since it provides a more reliable extrapolation. 

4.4. The Kennicutt-Schmidt law 

We used our data to calculate the Kennicutt-Schmidt law, i.e. 
the relation between the disk-averaged gas column density and 
the disk-averaged SFR per area, Zsfr- Fig. [TTh shows the rela- 
tion with the molecular gas column density, £h2, and Fig. [T7b 
with the total (molecular+atomic) gas column density, 2h2+hi- 
The surface densities were calculated by dividing the SFR, re- 
spectively the (extrapolated) molecular or total gas mass, by the 
galaxy surface ttDos 1^- ^ c ^ ear correlation exists with the molec- 
ular gas, but none with the total gas, showing a lack of corre- 
lation with the atomic gas column density. Our findings are in 
agre ement with previous results based on spatially resolved anal - 
ysis dWong & Blitdl2002tlBigiel et alJl2008t IVerlev et al.l l2010). 
showing that the SFR is strongly related to molecular gas only. 
The best-fit regression yields Ssfr k 2^ 89±()07 (Table [7J. This 
slope, close to unity , is again consi s tent w ith the results of 
IWong & Blitzl (120021) and lBigiel et all (120081) for a spatially re- 
solved analysis. 

5. Discussion 

We compare our data to that of other samples of "normal" galax- 
ies and also search for differences in samples of interacting 
galaxies. The samples of normal galaxies that we consider are 
made of not obviously interacting galaxies, which are however 
selected without a clear isolation criteria, and some of these sam- 
ples contain, e.g. several cluster galaxies. In all these compar- 
isons we have adjusted the molecular gas masses to our defi- 
nition, i.e. same conversion factor and no consideration of the 
helium fraction. 

5. 1 . Comparison to studies of normal galaxies 

5.1 .1 . Relation between molecular gas, FIR and blue 
luminosity 

The nonlinear relation betwe en Mh, and Lb has been found 
by other groups as well (see iPerea et aljfl997l and references 
therein), with similar slopes as found by us. Their study was 
based on a smaller number (n = 68) of galaxies selected with 
a less rigorous criterion with respect to the environment. They 
discussed the cause of the nonlinearity and conclude that the 
most likely reason is extinction affecting Lb and increasing with 
galaxy luminosity. They predicted that the relation between Mh, 
and luminosities at longer wavelengths should be more linear. 
This prediction is confirmed by the nearly linear relation found 
between Mh, and Lk in our analysis. 

We compared our mean value of Mh 2 /Lb for galaxies of type 
T = 3 - 5 to resu l ts for a large sample of normal galaxies, studied 
bv lBettoni et al.l ((2003). They searched the literature for galax- 
ies with data for their ISM properties, excluding galaxies with 
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Fig. 17. Panel a): The relation between the surface density of 
the molecular gas and the SFR per area for morphological types 
T = 3 — 5, calculated by dividing the molecular gas mass and the 
SFR derived from Lfir by the surface of the galaxies, n(D25/2) 2 . 
The full line shows the bisector best fit for this sample, and the 
dashed line the best fit to the sample of all morphological types. 
Panel b: The same relation with the total (atomic + molecular) 
gas surface density. The blue triangles show the galaxies with 
either an upper limit in HI or in CO. The gas surface density 
for these galaxies has a lower and an upper limit which for the 
galaxies in this figure are very close together (within the size of 
the triangle). 



The relation between Mh 2 and Lfir in nonstarburst galax- 
ies has been found to be close to linea r in other studies (e.g. 
iPerea et al]|1997HGao & Solo mon 2004), in agreement with our 
results. 



5.1 .2. Gas depletion time 

We derived in Sect. 4.2.4 the SFE (defined as the SFR/M H2 ). The 
SFE is directly related to the molecular gas depletion time, Tdep, 
by Tde P = SFE -1 . The gas depletion time for our sample has a 
mean value of log(Td ep ) = 8.9 yr for the entire redshift-limited 
sample and log(Tiep) = 9 yr for galaxies of type T = 3 - 5 
(see Table |6j, with a spread of values roughly ranging between 
log(rdep) = 8.5 yr and log(r dep ) = 9.5 yr (see Fig. [Till. 

This val u e can be compared to those found in recent surveys. 
Bigi el et alJ d20TTb derived a mean gas depletion time of 2.35 
Gyr from spatially resolved observations of 30 nearby galaxies 
in the HERACLES survey. This values includes a helium frac- 
tion of a factor 1.36, thus giving Td ep = 1.7 Gyr without helium. 
They furthermore showed that this value is consistent with a 
wide range of molecular gas depletions times from the literature, 
albeit with a large standard deviation of 0.23 dex. This value 
is only slightly higher than our value for isolated galaxies. The 
small difference could be due to the fact that our value is global, 
and thus might encompass some FIR emission not directly asso- 
ciate d to SF from the ou tskirts of the galaxies, whereas the value 
from Bigiel et al. ( 201 ll) is from a spatially resolved study. 

ISaintonge et al.1 d201 lb studied the molecular gas depletion 
time, Tdep, for a volume limited sample of 222 galaxies with 
0.025 < z < 0.05 observed in the COLD GASS survey. They 
found values for the gas depletion time in the range of roughly 
log(T d eD ) = 8.6 yr and 9.5 yr, with a mean value of Td ep = 1 Gyr 
ffor a IChabrierl (120031) IMF which gives a mass comparable to 
within ~ 10% to the Kroupa IMF). Both the range and the mean 
agree very well with our values. 

They furthermore found a good correlation of T dep with M„ 
and with the specific SFR, sSFR = SFR/M,. We tested both cor- 
relations with our data. We did not find a similarly good corre- 
lation with Li((which is in a good approximation proportional to 
the stellar mass M.) although our sample covered a similar range 
of stellar masses, up to several 10 M . Taking into account all 
morphological types, we found evidence for a weak trend (corre- 
lation coefficient r = 0.25) which completely disappeared when 
restricting the sample to T — 3 - 5 (r — 0.08). 

We could confirm with our data the existence of a correlation 
between Td ep and sSFR. This correlation is however expected be- 
cause both Tdep and the sSFR depend on the SFR, and Lk and 
Mh, show a good correlation. Thus, the respective ratios are ex- 
pected to correlate. 



a known pecularity (interacting, disturbed, galaxies with polar 
rings or counterrotation) and with active galactic nuclei. Their 
sample includes 177 galaxies of type T = 3 - 5 with CO(1-0) 
data (160 detections and 17 upper limits) with values of log(LB) 
between 9 and 11. They derive a mean value (adapted to our 
convention for the calculation of the Mh 2 and Lb) of MylJLq= 
-0.92 ± 0.04 for these galaxies. When restricting the range of Lb 
to log(L B ) = 10 - 10.6, the mean value is -0.82 + 0.05, about 0.2 
dex higher than the corresponding value for the AMIGA sample 
(see Tab. [8). showing that AMIGA galaxies have a lower molec- 
ular gas content. When comparing their values for Mh 2 /D~ a 
similar difference is found. 



5.1.3. The molecular-to-atomic gas ratio 

The value of MnJMm in galaxies as a function of morphological 
type has been cont roversial. In a large survey of spiral galaxies, 
Young et al. ( 1995) map ped a sample of about 3 00 galaxies with 
the FCRAO telescope. [Young & Knezekl (Tl989) studied the de- 
pendence of Mh, on Mhi for that sample. They derived a conti- 
nous decrease of M^/Mni from early to late-type galaxie s with 
mean values above 1 for early types. ICasoli et al.1 (U~998) stud- 
ied the molecular gas properties of a large (n = 528), heteroge- 
nous sample of galaxies, composed of data from the literature 
and their own observations. Their data consisted both of mapped 
galaxies and of objects where only the central position had been 
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observed. In agreement with lYoung et al.l (1 19951) . they found a 
decrease in Mh,/Mhi from early to late-type galaxies, but ob- 
tained much lower values for Mh,/Mhi, especially for early-type 
spirals. In Fig. [18] we compare our resul ts to these two st udies 
and furthermore to the THINGS s ample (|Lerov et aL 20081) and 
the No beyama survey dNishivama & NakafeOoilT Betton i et al.l 
(120031) (not included in the Fig. [18]) prov ide values fo r a sam ple 
of 427 normal galaxies in agreement with lCasoli et al.l d 1998b for 
T > and higher values (by 0.5-1 dex) for T < 0. 

The data from all studies show a decrease of the Mh,/Mhi gas 
mass towards late morphological types. A pronounced step in 
the ratio takes place at T ~ 6. The low values found for late-type 
galaxies co uld be due to two effects: ( i) late-type galaxies are 
richer in HI (lHavnes & Giovanellill984l) or (ii) they have a lower 
molecular gas content. A comparison of Mhi and Mh, to the 
blue luminosity or the optical diameter shows that both effects 
take place: The ratio Mhi/Lb (respectively Mhi/L) 2 ^) increases 
for types T > 6 by ~ 0.2 - 0.3 dex wher eas MuJLr (respectively 
M H ,/L>? 5 ) decreases by ~ 0.3 - 0.5 dex dYoung & Knezeklll989t 
ICasoli et afl 19981: IBettoni et al.ll2003l) . Our data show the same 
behaviour. The strong decrease of Mh,/Lb and M^JD 1 ^ could 
be due to a real decrease in the molecular gas mass or due to 
the fact that late-type galaxies tend to have lower metallicities 
so that we probably underestimate the true Mh, by using the 
Galactic conversion factor. 
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Although the general trends in the data sets are the same, 
with a pronounced and continous decrease of Mh 2 /Mhi from 
early to late-type galaxies, there are considerable differences be- 
tween the different samples (see Fig.lTSTi. The differences are up 
to an order of magnitude for elliptical, lenticular and early-type 
spiral galaxies (up to T ~ 3), and less (up to 0.4 dex) for later 
type spirals. There is in general a good agreement between the 
values found by Casoli et al. and ours, the match being very good 
for early-type galaxies (up to T= 0) whereas for later-type galax- 
ies our values he somewh at below. The results from the THINGS 
survey (|Lerov et al1l2008l) are also in agreement with the data of 
ICasoli et all and only sligh tly h igher than our values. The r esults 
bv lYoung & Knezekl d 19891) and lNishivama & Nakail d2001l) give 
much higher values than the other surv eys. This could be due to 
the sample selection. In the sample of lYoung & Knezekl |l989) 
FIR or optically bright galaxies were selected, so that a high CO 
emissi on can be expected. Furthermore, the I Young & Knezekl 
( 1989) survey contains objects from the Virgo Cluster, and clus- 
ter galaxies are known to b e HI defic ient on average. Cluster 
galaxies were excluded in Casol i et ail and are not present in 
our sample either. A selection effect can also account for the dif- 
ference between the samples of Nishiv ama & NakalllLerov et al.l 
and our data. The sa mples differ in the mean Lfir/Lb with the 
Nishivama & Nakai sample being the FIR brightes t with a mean 
log(LpiR/LB) = 0.10±0.07 whereas the lLerov et alJ and our sam- 
ple show lower values (log(LpiR/LB)= -0.51 ± 0.06, respec- 
tively log(LpiR/iB)= -0.49 + 0.03). Since galaxies with a higher 
logiLpiR/Lft) tend to have a higher Mh,/Mhi, the observed trend 
is understandable. 

A noticeable point is that the mean M H ,/M H i of the AMIGA 
sample represents for all morphological types the lower limit of 
all samples. For late-types, the AMIGA galaxies have the low- 
est Mh,/Mhi, even though the Mh, is extrapolated unlike other 
studies that are usually based on a central pointing only. This 
is in line with the low mole cular gas content th at we found in 
comparison to the sample of Bettoni et aD (12003b . 



5.2. Comparison to studies of interacting galaxies 

5.2.1. Is the molecular gas content enhanced in interacting 
galaxies? 

Several studies in the p ast have concluded that Mh, in interacting 
galax i es is enhanced (|Braine & Combesl Il993t ICasasola et al.l 
l2004t ICombes et al.lll994l) . based on a higher value of M H ,/L B 
or Mh 2 /£>25 compared to noninteracting galaxies. However, we 
found that the ratio Mh,/Lb increases with Lb and this trend has 
to be taken into acco unt when compari ng isolated and interact- 
ing samples. Indeed. IPerea et all d 1997b found no difference in 
the correlation between Mh, and Lb for samples of isolated, 
strongly and weakly perturbed galaxies. 

We use the sample of Casa sola et alJ d2004l) of interact- 
ing galaxies to search for a possible excess in Mh, with re- 
spect to AMIGA galaxies. Their sample includ es 153 galaxies 
with molecular gas data from different sources. ICasasola et al.l 
found that the mean ratio of Mh,/Lb for spiral galaxies was, 
depending on the morphological type, between about 0.2 and 
1.0 dex hi gher than for a sam ple of 427 noninteracting galax- 
ies from IBettoni et al.1 d2003l) and concluded that Mh, was 
enhanced in interacting galaxies. However, since their sam- 
ple of interacting galaxies is on average 0.5 mag brighter, a 
higher Mh,/Lb is already expected due to the higher luminos- 
ity. We furthermore inclu de the sampl e s of s trongly and weakly 
perturbed galaxies from IPerea et al. I (fl997h in this test. The 
weakly perturbe d sample has 43 g alaxie s and includes classes 
1, 2, and 3 of ISolomon & S age (1988) a nd cla s s 2 o bjects 
from the luminous IRAS sample of ISanders et al.l d!991l) . The 
strongly perturb ed sample has 35 galaxi es and includes inter- 
action cla ss 4 of ISolomon & Sagd (ll988). interaction classes 3 
and 4 of Sanders et al. (1991) and closely interacting pairs from 
ICombeset al.ldl994l) . 

We searched for a possible enhancement of Mh, in compar- 
ison to Lb and to Lk in these samples, both by applying the de- 
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ficiency parameter and by comparing the ratios. The mean val- 
ues are listed in Table [8] We note that the interacting samples 
are more luminous (by about 0.5 dex) in both Lb and Lk than 
the AMIGA sample so that we have to extrapolate the relations 
found for the AMIGA galaxies to higher luminosities. For the 
calculation of the deficiency parameter we therefore use the bi- 
sector fit as the best fitting method for extrapolations. We use 
the fits for the AMIGA T — 3 - 5 subsample to compare to 
the weakly and strongly perturbed samples WPER and SPER, 
since spiral galaxies are the b etter compariso n for these actively 
star-forming objects. For the lCasasola et al.l sample information 
about the morphological types are available and we are able to do 
the analysis both for the entire sample and for T = 3-5, using the 
values for the corresponding morphological types in AMIGA. 
When comparing the ratio Mh,/Lb we restrict the galaxies that 
we consider to the same range of Lb in order to avoid effects 
caused by the nonlinearity in the correlation. For the roughly 
linear Mh,/Lk ratio this is not necessary (we checked that no 
difference in Mh,/Lk was present for low and high Lk in the 
AMIGA sample). 

Fig[19] shows Mh, vs. Lb for the three samples compared to 
the bisector fit of the AMIGA T — 3 - 5 sample. An excess 
of Mh, is visible for both the Casasola T — 3 - 5 and SPER 
sample. This is confir med by the mea n value of the deficiency 
(Table H), both for the ICasasola et al.l T = 3 - 5 and the SPER 
sample a clear excess (i.e. a negative value of the deficiency) of 
Mh, is found. This excess in the SPER sample was not found 
in Per ea et al.l dl997h due to the different regression line they 
used. For the WPER sample, on the other hand, no deficiency nor 
excess is found. When comparing the ratio Mh,/Lb we confirm 
the higher values for both the Casasola and the SPER sample, 
with mean values of about 0.3 dex higher than for the AMIGA 
T — 3 - 5 sample. The value of Mh,/Lb for the WPER sample 
is compatible with AMI GA within the errors. The entire (i.e. 
all morphological types) Casasol a" et al.l sample has a value of 
Mh,/Lb compatible with AMIGA, showing the importance of 
taking the morphological type into account. 

When doing a similar analysis based on Lk (see Fig. [20b . 
we find again no indications for an enhancement in Mh, for 
the WPER sample nor for the total Casasola sample, neither 
from the deficiency nor from the ratio. There is a clear excess 
in Mh, for galaxies of type T — 3 - 5 from the Casasola sam- 
ple, quantified both in the ratio, which is about 0.3 dex higher 
for T — 3 - 5 AMIGA galaxies, and in the deficiency parame- 
ter. Also the SPER sample shows a higher Mh,/Lk ratio than the 
AMIGA sample, however, no indication of an Mh, excess from 
the deficiency parameter. The discrepancy between the two indi- 
cators (mean and deficiency) is probably due to the small sample 
size and the small amount of the molecular gas excess. We sug- 
gest to rely in this case on Mh 2 /Lk as the more robust indicator. 

In summary, by comparing the molecular gas mass to Lk 
and Lb and using both the deficiency and the ratios Mh,/Lk 
and Mh,/Lb we found the clearest evidence for an enhancement 
of Mh, in the Casasola sample of interacting galaxies of type 
T = 3 — 5, which is the largest comparison sample in our study. 
Evidence for an enhancement was also found for the Perea sam- 
ple of strongly perturbed interacting galaxies (SPER). No ev- 
idence for any enhancement was found for the Perea sample 
of weakly perturbed interacting sample (WPER). Based on the 
present data, it is not entirely clear where the differences between 
these samples come from. A possible reason could be a lower 
degree of interaction in the WPER sample. We would like to 
point out the importance of matching the luminosity range when 
comparing parameters with a nonlinear correlation, as Mh,/Lb, 




9 10 11 

log(L B ) [L ] 

Fig. 19. The molecular gas mass vs. Lb for different sample of 
in teracting galaxi e s: Gal axies of type T — 3 - 5 from the sample 
of ICasasola et al.l (|2004) and a sample of weakly (WPER ) and 
strongly (SPER) perturbed galaxies from lPerea et al.l(ll997l) . We 
adapted both the molecular gas masses and Lb to our definition. 
The full line is the regression fit obtained for the AMIGA T = 
3-5 sample from Table [7] 



and the importance of comparing the same morphological types, 
since there are generally large difference between early type and 
spiral galaxies, in particular for Mh,/Lk. 

What influence has the fact that the molecular gas masses in 
these samples are not aperture corrected as in the AMIGA sam- 
ple? We do not expect this effect to be very important in inter- 
acting galaxies where the SF and thus the molecular gas usually 
tend to be more concentrated to the central regions. However, 
we can make an e stimate of the importance of this effect for 
the ICasasola et al.l sample. The median angular diameter of the 
galaxies of type T — 3 - 5 in this sample is 90". If we assume 
that the galaxies were observed only at the central position with 
a beam size of 50" (the typical beam size of the radio telescopes 
used in surveys) the predicted aperture correction according to 
our prescription is between a factor of 1 .2 and 1 .4 (for edge-on 
and face-on galaxies, respectively). Thus, in this case the total 
molecular gas content would be 0.08-0.15 dex higher. Thus, if 
the molecular gas in the interacting galaxies of the Casasola sam- 
ple are distributed in a similar way as in isolated galaxies, this 
difference would strengthen the finding that Mh, is enhanced in 
this sample of interacting galaxies. It also shows the importance 
of mapping the molecular gas in galaxies in order to be able to 
compare different samples in a relyable way. 

5.2.2. Is Lfir and the SFE enhanced in interacting galaxies? 

We use the same samples as in the previous subsection to look 
for a possible enhancement of Lfir in interacting galaxies. In 
Fig. [2T] Lfir is compared to the blue luminosity. We include the 
best-fit regre ssion line found for th e AMIGA T — 3 - 5 sample 
presented in lLisenfeld et al.l (120071) (Table 6, log(L F iR) = 1.35x 
log(LB)-3.98). We see a clear excess of Lfir with respect to this 
regression line for the three samples, the strongest for the sample 
of strongly perturbed galaxies (SPER). The excess can be quan- 
tified by calculating the Lfir deficiency, defined in an anologous 
way as the Mh, deficiency. The values are listed in Table [8] and 
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AMIGA Casasola et al. (2004) Perea et al. (1997) 





all types 


T = 3-5 


all types 


T = 3-5 


WPER 


SPER 


< def(M H , )> (L B , bisec.) 


0.07±0.04 


0.07±0.04 


0.13+0.07 


-0.31±0.07 


-0.03+0.07 


-0.25+0.08 


ii/n uf 


173/79 


88/21 


153/35 


68/9 


43/0 


35/0 


<log(M H2 /L B )> (for 10 10 <L B > 10 10 ' 6 ) 


-1.10±0.06 


-1.04±0.05 


-1.12±0.08 


-0.71±0.08 


-0.93+0.06 


-0.71+0.08 




50/11 


40/5 


85/13 


39/6 


24/0 


19/0 


< def(M H;! )> (L K , bisec.) 


-0.14±0.19 


0.05±0.04 


0.03±0.08 


-0.16±0.06 


0.13+0.07 


-0.02+0.05 


<log(M H2 /L K )> 


-1.87±0.05 


-1.66±0.06 


-1.70±0.08 


-1.37±0.07 


-1.65+0.07 


-1.47+0.05 


»/«up 


173/79 


88/21 


132/29 


60/6 


43/0 


32/0 


<def(L P1R )> 


-0.02±0.04 


-0.05±0.04 


-0.39±0.02 


-0.44±0.03 


-0.48+0.08 


-1.07+0.09 


n/n up 


172/75 


88/20 


628/0 


340/0 


43/0 


35/0 


<log(L F1R /M H2 ) > 


0.72±0.03 


0.63±0.03 


0.89±0.05 


0.75±0.05 


0.98+0.07 


1.33+0.05 


n/n up 


97/22 


68/10 


105/15 


51/4 


43/0 


35/0 
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Fig. 20. The molecular gas mass (adapted to our conversion fac- 
tor) vs. Lk for the same samples as in Fig. [19]. The full line is the 
regression fit obtained for the AMIGA T — 3 - 5 sample from 
Table |7] 



confirm the visual impressions of an excess in Lfir of about an 
order of magintude for the SPER sample. 

The large excess in Lfir together with a smaller (or no) ex- 
cess in Mh 2 results in a higher value of SFE oc L fir /Mh, com- 
pared to the AMIGA sample for the strongly and weakly per- 
turbed samples. Comparing these values to the results for the 
AMIGA sample we find an increase of about a factor 5 for the 
strongly perturbed sample and of about 2 for the weakly per- 
turbed sample. The SFE in the Casasola sample is similar to the 
AMIGA sample. 

For higher infrared lumi nosities, the v a lue of Lfir/Mh, is 
known to increase strongly. ISanders et alJ (Il99ll) studied the 
molecular gas content in luminous infrared galaxies, ranging 
from Lir = 10 10 L to several 10 12 L , and showed that the 
ratio Lir/Mh 2 increases strongly with IR luminosity and with 
the degree of interaction. Whereas for isolated, low-luminosity 
(Lir< 10 11 L Q ) galaxies they found values of Lir/Mh 2 similar to 
ours, Lir/Mh, increases by a factor of about 10 for galaxies with 
Lir~ 10 12 L Q which are mostly advanced mer gers. Similarly 
high v alues of Lfir/Mh, ~ 50 have been found bv lSolomon et alJ 
(1997) for ULIRGs. 
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A Perea-SPER 
• Perea-WPER 
x Casasola 




Fig. 21. The FIR luminosity v s. blue luminosity fo r different 
samples of interacting galaxies: ICasasola et al.l ((2004) (only spi- 
ral galaxies with type T — 3 - 5 are shown, in order not to over- 
load the plot), and a sample o f weakly (W PER) and strongly 
(SPER) perturbed galaxies from lPerea et al.l dl997l) . The full line 
is the regression fit obtaine d for the AMIGA T = 3 - 5 sample 
from Lisenfeld et al.l (120071) . 



6. Conclusions and Summary 

We presented molecular gas masses, based on CO observations, 
for a sample of 273 isolated galaxies and we performed a statisti- 
cal analysis for a redshift-limited sample of 173 isolated galaxies 
with recession velocities between 1500 and 5000 km s _1 . The 
observations covered in most cases only the central position of 
the galaxies. In order to correct for the missing molecular gas 
mass outside the observed area, we derived and applied an aper- 
ture correction assuming an exponentially decreasing CO disk. 
We obtained the following results: 

1 . We compared the molecular gas mass to different parameters 
(Lb, Lk, and Lfir) in order to characterize the relations 
followed by isolated galaxies and establish a baseline that 
can be used to find possible deviations for interacting galax- 
ies. We concentrated our analysis on the morphological types 
T = 3 - 5 (Sb-Sc) which represents the bulk of our sample 
(51% of the galaxies are of these types) where also the de- 
tection rate of CO (74%) was highest. 
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2. We found good correlations with these parameters, roughly 
linear in the case of Lk and Lfir and nonlinear for Lb and 
Dj 5 . The tightest correlations are with Lfir and, for spiral 
galaxies with T = 3 — 5, with Lk, whereas the correlation 
with is the poorest. Due to the nonlinearity of the cor- 
relation, the ratio MnJ Lb changes with Lb which has to 
be taken into account when comparing it to other samples. 
We describe a deficiency parameter, defined in analogy to 
the deficiency parameter for the atomic gas as the difference 
between the logarithm of the expected molecular gas mass 
and the logarithm of the observed molecular gas mass. The 
expected molecular gas mass can be calculated from any of 
the parameters studied by us (Lb, D25, Lk, Lfir) using the 
correlation coefficients listed in Table|7] 

3. We applied these relations and the res ulting expressions 
to three samples f rom the literature dPerea et alJ Il997t 
ICasasola et alj|2004h . For the sample of Casasola et al. and 
the sample of strongly interacting galaxies of Perea et al., we 
found clear evidence for an enhancement of Mh, in compar- 
ison to Lb and Lk, while for a sample of weakly interacting 
galaxies from Perea et al. no difference with respect to the 
AMIGA sample was found. A possible reason for this differ- 
ence could be a higher degree of interaction in the first two 
samples. 

4. We derived a mean molecular gas depletion time, Td ep (de- 
fined as Mh,/SFR), of log(Tdep) = 9.0 yr for spiral galaxies 
(T = 3-5) and a slightly lower value of log(Td ep ) = 8.9 yr for 
all morphological types (both values for a Kroupa IMF), in 
reasonable agreement with other s tudies of nearby galaxies 
dBigiel et alj201 UlSaintonge et alj|201 ll) . 

5. No good correlation was found between Mh 2 and Mm. The 
ratio between the molecular and the atomic gas mass de- 
creases significantly from early to late-type galaxies, with 
a difference of up to a factor of 10. The ratio MuJMm of the 
AMIGA galaxies is well below 1 for all morphological types, 
with a mean value of log(MH 2 /Mni) = -0.72 for galaxies of 
type T = 3-5. We compared our values to those of other 
noninteracting samples and found that the AMIGA galaxies 
had the lowest values for all spiral galaxies. 

6. We used our data to compare the disk averaged surface 
densities of the molecular, Sh, , and molecular+atomic gas, 
Xh 2 +hi, to those of the SFR, Ssfr- We found a good correla- 
tion between the logarithms of Sh 2 and Xsfr, with a slope 
close to 1. No correlation with of Ssfr with £h 2 +hi was 
found. 
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Appendix A: CO Spectra 



Fig. lA.ll shows the CO(1-0) profiles of the detections and tenta- 
tive detections observed by us at the IRAM 30m telescope and 
Fig. lA.2l those observed at the FCRAO 14m telescope. 
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Fig.A.l. CO(l-O) spectra for the galaxies detected with the IRAM 30m telescope. The x-axis represents the recession velocity 
in km s _1 and the y-axis the main beam brightness temperature T m \, in K. The spectral resolution is 20.8 km s _1 in most cases 
except for some individual galaxies for which a higher or lower resolution was required to clearly show the line. The full (red) line 
segment shows the line width of the CO line adopted for the determination of the velocity integrated intensity. The dashed (black) 
line segment is the HI line width at 20% peak level, W2o- An asterisk next to the name indicates a marginal detection. 
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Fig. A.l. (continued) 
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U. Lisenfeld et al.: The AMIGA sample of isolated galaxies. X 




Fig. A.l. (continued) 
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Fig.A.2. CO(l-O) spectra of the galaxies detected with the FCRAO telescope. The spectral resolution is 13.1 or 26 km s _1 , de- 
pending on the spectral shape and rms noise. The x-axis represents the recession velocity in km s _1 and the y-axis the intensity in 
K in the T m b scale. The full (red) line segment shows the line width of the CO line adopted for the determination of the velocity 
integrated intensity. The dashed (black) line segment is the HI line width at 20% peak level, W2o- An asterisk next to the name 
indicates a marginal detection. 
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Fig. A.2. (continued) 
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